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Abstract 

We have investigated the finite temperature elastic properties of AIRE (RE=Y, Tb, Pr, Nd, Dy) with B2-type 
structures from first principles. The phonon free energy and thermal expansion is obtained from the quasi- 
harmonic approach based on density-functional perturbation theory. The static volume-dependent elastic 
constants are obtained from energy-strain functions by using the first-principles total-energy method. The 
comparison between our predicted results and the ultrasonic experimental data for a benchmark material Al 
provides excellent agreements. At T = OK, our calculated values of lattice equilibrium volume and elastic 
moduli of our calculated AIRE (RE=Y, Tb, Pr, Nd, Dy) intermetallics agree well with the previous theoret- 
ical results. The temperature dependent elastic constants exhibit a normal behavior with temperature, i.e., 
decrease and approach linearity at higher temperature and zero slope around zero temperature. Furthermore, 
the anisotropy ratio and sound velocities as a function of temperature has also been discussed. 
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1 Introduction 

Aluminium alloys are widely used in the automotive, aircraft and aerospace industries, and interest of their applica- 
tions is rapidly increasing. Addition of rare earth (RE) elements to Al-based metal alloys have received great attention 
since it can result in increasing their wear resistance 131, mechanical properties |j2l, electrochemical behavior [3| and 
thermal stability |4|,etc. The development and design of new Aluminum alloys need more fundamental physical data 
of these alloys. Recently, some investigations focusing on AIRE, AI2RE, and AI3RE systems intensely are performed 
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from first-principle calculations based on density functional theory (DFT) lH] |7] |6l [8l |9l [lOl [TT] [T2| . The mechanical and 
thermodynamical properties of B2-A1RE intermetallics have been investigated from using the projector augmented wave 
(PAW) method within generalized gradient approximation (GGA) by Tao and Ouyang et al EH, while Srivastava et al 
lis] also investigated the electronic and thermal properties using the tight binding linear muffin tin orbital (TB-LMTO). 
These results are very important to the novel material design and further scientific and technical investigations. 

Elastic properties of a solid are very important because some physical properties such as the bulk modulus, shear mod- 
ulus. Youngs modulus, anisotropy, and the velocity of sound can be determined from its elastic constants IIT3l . Further- 
more, the finite temperature elastic constants are critical for predicting and understanding the mechanical strength, stabil- 
ity, and phase transitions of a material lfT4l . Though modem electronic-structure methods based on DFT can compute the 
zero-temperature pressure dependence elastic moduli from first principles accurately ifTSl . treating the corresponding tem- 
perature dependence of elastic constants is still a formidable challenge lfT6l . However, several groups have also attempted 
to investigate the temperature-dependent elastic properties by using the first-principles DFT methods 1161 fTTl [TSl [191 l20l . 
More recently, we have investigated the temperature-dependent elastic properties of two ductile rare-earth intermetallics 
YAg and YCu, and MgRE intermetallics from using the first principles approach successfully 11211 1221. So far as we know, 
the temperature-dependent elastic moduli of the novel B2-tpye structures AIRE intermetallics have not been reported in 
the literature. We investigate the finite temperature elastic properties for AIRE (RE=Y, Tb, Pr, Nd, Dy) using a first- 
principles method |fT9l l20l in present study. The elastic constants as a function of temperature, the anisotropy ratio, and 
the temperature-dependent sound velocity have been discussed. 

2 Theoretical methods 

The isothermal elastic constants are derived from the second-order strain derivative of the Helmhotz free energy ||23]| 



where Q is equilibrium volume at temperature T, s denotes the linear strain tensor, and the Helmholtz free energy F can 
be expressed as the function of volume Q. and temperature T, 



where is the zero-temperature total energy of the electronic ground state, F^i represents the electron-thermal con- 
tribution from finite temperature electrons, and Fph denotes the phonon free energy arising from the lattice vibrations. 
For the electron-thermal free energy, an acceptable approximation to calculate it is to use the Mermin statistics ll24l |25]| 
through Fei - Ef.\ - TSei. Within the mean-filed approximation, the thermal electronic energy Eel and entropy 5ei 
can be give by E,, = / e«(6)/(6)rfe - f " e«(e)^/e and 5ei = -^b / «(e)[/(e) ln/(e) + (1 - /(e))ln(l - /(e))]t/e, re- 




■ijOSki T,s=Q 



(1) 



T) = EoiD.) + Fei(n, r) + Fph(Q, T), 



(2) 
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spectively, where fcg is Boltzmann constant, ef is Fermi level, /(e) is the Fermi-Dirac distribution, and n(e) is the den- 
sity of states. The phonon free energy can be obtained within the framework of quasiharmonic approach (QHA) by 
Fph = Yjk \\^^k + ksT In (1 - e-'''^«/*B^)j^ where represents an individual phonon frequency. 

For the B2-type AIRE intermetallics with cubic crystal, there are only three independent constants C[^|, Cjj, and CJ4 
(in Voigt notation). A more convenient set for computations are Cj^ and two linear combinations Bj and fij- The bulk 
modulus Bj - {C\y + IC^^)/^ is the resistance to deformation by a uniform hydrostatic pressure and can be obtained 
from the the equation of state (EOS) ||T41 . using the Vinet EOS 1261 in the current calculation. The shear modulus 
fiT = (C^j - C[2)/2 is the resistance to shear deformation across the (1 10) plane in the [110]. 

We apply a volume conserving tetragonal strain to calculate shear modulus fir, 

e(S) ^ (6,6,(1 +Sr^ -1,0, 0, 0), (3) 

and the corresponding free energy is 

F(n, 6) = F(n, 0) + 6^1x^6^ + 0(6^), (4) 
where F(Q, 0) is the free energy of the unstrained state. Cjj is obtained from a volume-conserving monoclinic stain, 

b(6)^ (0,0,6^ /(I -6^), 0,0,26) (5) 

which leads to the free energy change 

F(Q, (5) = F(Q., 0) + 2Cl^Q.6^ + 0(6''). (6) 

The elastic constants are calculated from the fitted quadratic coefficients of the strain free energy. Then, the elastic 
constants C[^| and Cjj are obtained from jij and Bj. 

The electronic energy, including the static energy at ground state and thermal electronic contribution to the free energy, 
are computed by using the density functional theory (DFT). For the first-principles DFT calculations, we employed the 
plane-wave basis projector augmented wave (PAW) method Il27l l28l in the framework of of the generalized gradient 
approximation (GGA) in the Perdew-Burke-Ernzerhof form |l29l|30l as implemented in the VASP code CT [32l l33l [34l . 
Since high accuracy is needed to evaluate the elastic constants, the convergence of strain energies with respect to the 
Brillouin zone integration was carefully checked by repeating the calculations for 21 x21 x21 and 25 x 25 x 25 Monkhorst- 
Pack ||35l meshes at equilibrium volume at T=OK, and we found that fluctuation both for CJ4 and jij is lower than 0.5GPa. 
Hence, we used 21 x21 x21 in the full Brillouin zone giving 726 irreducible k-points. In addition, we used a high plane- 
wave energy cutoff of 600eV which is sufficient to calculate the elastic moduli accurately. 

To calculate the phonon free energy in Eq. Q, we have carried out the density-functional perturbation theory (DFPT) 
within supercell approach. Real-space force constants of supercells are calculated in the DFPT implemented in the VASP 
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code ll34ll . and phonon frequencies are calculated from the force constants using the PHONOPY package ll36l [37l [38l . 
The chosen supercell size strongly influences on the thermal properties. We give Fpf, at 300K and lOOOK with 2x2x2, 
3x3x3, and 4x4x4 supercells for AIRE (RE=Y, Tb, Pr, Nd, Dy) in Table [T] We found that the phonon free energies 
between 3x3x3 and 4 x 4 x 4 are less that 0.5kJ/mol, so the adequate supercell size consisting of 3 x 3 x 3 unit cells 
containing 54 atoms is used. The total Helmholtz free energy including the electronic energy and phonon free energy are 
calculated at 13 volume points. The equilibrium volume Q.q{T) at temperature T and the bulk modulus Bj are obtained by 
minimizing Helmholtz free energy Eq. (j2]) with respect to O from fitting the integral form of the Vinet equation of state 

(EOS) ma. 

The temperature-dependent elastic constants are obtained in the following three steps |fT9l l20l . The first step is 
calculating the the equilibrium volume D.o and thermal expansion as a function of temperature T by using the first- 
principles quasiharmonic approximation (QHA), while isothermal bulk modulus Bt is also obtained from fitting to the 
Vinet equation of state |26|. The detailed procedure of calculating the thermal expansion coefficient, isothermal bulk 
modulus Bt were shown in our recent work ||39l . In the second step, we predict the static elastic constants as a function 
of volume Cjj{Q.()) by using the energy-strain relation based on Eqs.(|3jj6jl. In the third step, the calculated elastic constants 
from the second step at the volume Qo(r) are approximated as those at finite temperatures, i.e., CfXT) = Cf:[T{Q.o)]. 

U 'J 

3 Results and discussion 

In Table |2] we list the equilibrium volume Q^, elastic constants C,j atT- OK for AIRE (RE=Y, Tb, Pr, Nd, Dy) in 
comparison with the previous calculated results obtained by Tao and Ouyang, et al.f7|. One can see that the values in our 
calculations are in excellent agreement with the previous DPT data |7 1. It is noticeable that the requirement of mechanical 
stability for all calculated AIRE, namely, Cn - C12 > 0, Cn > 0, and C44 > 1401 , is satisfied. In addition, the presented 
results of Al Y and AlPr are also agree well with theoretical results obtained from the tight binding linear muffin tin orbital 
(TB-LMTO) by Srivastava et al 15). 

The temperature-dependent elastic constants have not been reported theoretically and experimentally for all AIRE 
intermetallics, so we choose a benchmark material Al, which has the face-center-cubic (FCC) structure, accompanied 
by available experimental data to get a useful test of the accuracy of the method and the precision of our calculations. 
Furthermore, it is note that Eq. (|TJ gives the isothermal elastic constants Cjj 11411 . but most of the experimental data of 
elastic constants are usually measured as isentropic elastic constants. In order to compare with experimental data, the 
isothermal elastic moduli Cjj must also be transformed to the isentropic elastic moduli by the following relation ll42l 

Cu-Cl+^J]"^a,ClC], (7) 

kk' 

where Cq is the specific heat at constant volume and is the linear thermal expansion tensor. For the cubic crystal, Eq. 
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(|7]l simplifies to lfT6l 

with a is the volume thermal expansion coefficient, which is defined as = (l/Q.)dQ./dT. For the temperature depen- 
dence of the isentropic elastic moduli for Al and AIRE (RE=Y, Tb, Pr, Nd, Dy) intermetallics, we present in Figure [T] 
our calculated results. In Figure [TJ a), we also show our current date of Al together with data taken from ultrasonic mea- 
surements B3ll44ll . The overall agreement of the present results with experiment is very good at all temperatures, even 
near melt. Our results also show no discrepancy with previous calculated values obtained from using the first-principles 
approach by Wang et. al. |fT9ll . For all calculated AIRE intermetallics, we find the isentropic elastic constants C^j, ^^'^ 

decrease with increasing temperature, since thermal expansion may soften the elastic moduli. We also find that the 
trend of C^j , Cjj and approach linearity at higher temperature and zero slop around zero temperature. The calculated 
curves from Figure [T| indie ate that the elastic constants of AIRE follow a normal behavior with temperature lfT9l [T6l . 

The anisotropy ratio is defined as 

A = ^, (9) 

here we don't distinguish the isothermal and isentropic index for C44 and fi since they have same values. At T=OK, the 
anisotropy ratio has the values of 4.32, 4.40, 3.49, 3.71, and 4.48 for AIY, AlTb, AlPr, AlNd, and AlDy, respectively. The 
elastic anisotropy are estimated and it is shown that all the B2-A1RE compounds exhibit a substantial elastic anisotropy 
at K. In Figure [2] we present the anisotropy ratio as a function of temperature. We can see that the anisotropy ratios of 
AlPr, AlNd, and AlDy decrease linearly with increasing temperature, but those of AIY and AlTb decrease initially and 
increase after certain temperature, i.e., there are the minimum points in the temperature curves. This phenomenon can be 
understood from the temperature-dependent elastic constants in Figure[T] For AlPr [Figure [TJd)], AlNd [Figure [TJe)], and 
AUDy [Figure [ijf)], we find that the curves of C'^j, Cjj and Cf^ are nearly parallel, i.e., almost decrease to same extent in 
the whole temperature range. But for AIY [Figure[I|b)] and AlTb [Figure [TJc)], C^j decrease to slightly larger than 
and C44, and the differences among the elastic moduli become small. 

In addition, it is convenient to discuss the sound velocities as a function of temperature after obtaining the temperature- 
dependent elastic constants. From the ChristofFel equation ll45l . sound velocities are related to the elastic constants as 

(Cfji^,ninj-pv^6ij)ui^0, (10) 

where n is the propagation direction, u is the polarization vector, p is the density, and v is the velocity. For B2 intermetallics 
with cubic structure, the longitudinal mode along the [110] wave propagation direction is [I4i 

pvi = (C^i + Cfj + 2C^4)/2, (11) 



and two transverse modes are 

pv^i = Cl (12) 

and 

P^i = (Cn - Ci2)/2 (13) 

polarized along the [001] and [110] directions, respectively. For example, we show the sound velocities of AIY propagat- 
ing along the [1 10] direction as a function of temperature in Figure[3] Our results show that the sound velocities decrease 
with increasing temperature. However, the sound velocities vary slowly in comparison with the elastic constants since the 
values of sound velocities are proportional to the subduplicate elastic moduli. We also list the sound velocities at the room 
temperature {T = 3Q0K) in Table|3] and those data at finite temperature are very valuable for investigating the applications 
of AIRE intermetaUics. Owing to the strong anisotropy, two transverse modes along [110] direction split obviously. In 
the high-symmetry direction F - M, the order of V£ > vyi > vj-2 can also be obtained from the phonon dispersion curves 
directly. 

4 Conclusions 

In conclusions, the finite temperature elastic properties of AIRE (RE=Y, Tb, Pr, Nd, Dy) intermetaUics with B2-type 
structures from first principles method, in which the elastic constants as a function of temperature is obtained from the 
combination of the static elastic constants computed from energy-strain relations by the first principles total-energy cal- 
culations and the thermal expansion obtained from the first-principles phonon theory within the quasiharmonic approach 
based on density-functional perturbation theory. At OK, the equilibrium volumes, elastic constants for AIRE (RE=Y, Tb, 
Pr, Nd, Dy) in comparison with the previous calculated results are in excellent agreement with the previous DFT data. 
The overall agreement of the present results of the benchmark material Al with the ultrasonic experimental data is very 
good at all temperatures. The temperature dependent elastic constants of AIRE (RE=Y, Tb, Pr, Nd, Dy) intermetaUics 
exhibit a normal behavior with temperature, i.e., decrease and approach linearity at higher temperature and zero slope 
around zero temperature. The anisotropy ratios of AlPr, AlNd, and AlDy decrease linearly with increasing temperature, 
but those of AIY and AlTb decrease initially and increase after certain temperature, i.e., there are the minimum points in 
the temperature curves. Finally, we have discussed the sound velocities as a function of temperature, which shows de- 
creases with increasing temperature. However, the sound velocities vary slowly in comparison with the elastic constants 
since the values of sound velocities are proportional to the subduplicate elastic moduli. 
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Table 1: The phonon free energies at 300K and lOOOK with respect to supercell size. The unit of the energy is 



kJ/mol. 



Supercell size 


AIY 


AlTb 


AlPr 


AlNd 


AlDy 




300K 










2x2x2 


-5.032 


-7.393 


-7.684 


-7.605 


-7.584 


3x3x3 


-4.333 


-6.566 


-6.796 


-6.694 


-6.576 


4x4x4 


-4.384 


-6.546 


-6.787 


-6.682 


-6.556 




lOOOK 










2x2x2 


-78.266 


-85.919 


-86.789 


-86.559 


-85.598 


3x3x3 


-75.956 


-83.228 


-84.348 


-84.243 


-83.274 


4x4x4 


-76.137 


-83.109 


-84.124 


-84.109 


-83.115 



Table 2: The calculated equilibrium volume Q, elastic constants C,y, and bulk modulus B for AIRE (RE^Y, Dy, 
Pr, Sc, Tb) intermetallics at 7 = OK compared to previous first principles results. 





AIY 


AlTb 


AlPr 


AlNd 


AlDy 


a(A3) 


46.889", 46.850'' 


47.202", 47.202" 


53.157", 53.115'' 


51.853", 51.812'' 


46.539", 46.501'^ 


Cii (GPa) 


81.26", 77.93'' 


80.98", 77. 13'^ 


72.61", 66.83'' 


74.25", 69.00'' 


81.40", 78.28'' 


Ci2 (GPa) 


54.58", 56.58'' 


55.17", 56.92'' 


46.82", 49.26'' 


48.71", 51.09'' 


56.11", 57.73'' 


C44 (GPa) 


62.79", 61.91'' 


62.30", 61.17'' 


47.30", 46.67'' 


50.24", 49.83'' 


64.00", 63.06'' 


B(GPa) 


63.47", 63.70'^ 


63.73", 63.66'^ 


55.42", 55.11'' 


57.22", 57.06'' 


64.54", 64.58'' 



"This work 
''Referenced 



Table 3: The sound velocities propagating along the [110] direction including one longitudinal mode and two 
transverse modes at T=300K are listed. The unit is km/s. 





AIY 


AlTb 


AlPr 


AlNd 


AlDy 


vl 


5.573 


4.430 


4.487 


4.488 


4.397 


Vti 


3.894 


3.077 


3.027 


3.046 


3.067 


Vj2 


1.894 


1.488 


1.658 


1.612 


1.479 
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(f) AIDy 
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Figure 1 : The calculated isentropic elastic constants for (a) Al, (b) AIY, (c) AITb, (d) AlPr, (e) AINd, and (f) 
AIDy. The solid, dashed-dotted and dashed curves represent the values of C^j, and C^^, respectively. For 
Al, the empty symbols denote the corresponding experimental data of Kamm et al Il44l . and the full symbols 
indicate the experimental values from Gerlich et al [43j. 
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Temperature [K] 

Figure 2: The anisotropy ratio of elastic constants of AIRE (RE=Y, Tb, Pr, Nd, Dy) as a function of temperature. 
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200 400 600 800 1000 1200 
Temperature [K] 

Figure 3: The temperature-dependent sound velocity along the [110] propagation direction for AIY. The solid 
lines denote the longitudinal mode. The dashed and dashed-dotted curves represent two transverse sound 
velocities polarized along [001] and [lIO] directions, respectively. 
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